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ABSTRACT 

Coronal magnetic field may be characterized by how its field lines intercon- 
nect regions of opposing photospheric flux - its connectivity. Connectivity can 
be quantified as the net flux connecting pairs of opposing regions, once such re- 
gions are identified. One existing algorithm will partition a typical active region 
into a number of unipolar regions ranging from a few dozen to a few hundred, 
depending on algorithmic parameters. This work explores how the properties 
of the partitions depend on some algorithmic parameters, and how connectivity 
depends on the coarseness of partitioning for one particular active region magne- 
togram. We find the number of connections among them scales with the number 
of regions even as the number of possible connections scales with its square. 

There are several methods of generating a coronal field, even a potential field. 
The field may be computed inside conducting boundaries or over an infinite half- 
space. For computation of connectivity, the unipolar regions may be replaced by 
point sources or the exact magnetogram may be used as a lower boundary con- 
dition. Our investigation shows that the connectivities from these various fields 
differ only slightly - no more than 15%. The greatest difference is between fields 
within conducting walls and those in the half-space. Their connectivities grow 
more different as finer partitioning creates more source regions. This also gives a 
quantitative means of establishing how far away conducting boundaries must be 
placed in order not to significantly affect the extrapolation. For identical outer 
boundaries, the use of point sources instead of the exact magnetogram makes 
a smaller difference in connectivity: typically 6% independent of the number of 
source regions. 

Subject headings: MHD — Sun: corona — Sun: magnetic fields 
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Introduction 



According to prevailing understanding, coronal activity on the Sun involves energy 
stored in its magnetic field. The coronal field is stressed as the photospheric regions to 
which it is anchored slowly evolve. In order to build a quantitative model based on this 
insight it is essential to quantify how the coronal field links these photospheric regions - its 
connectivity. 

Several previous studies have focused on local properties of magnetic connec tivity as 
characterized by the point-wise mapping of positive footpoints to n egative footpoints (ISeehafei 
19861 : E3[l983; llnveraritv fc Titovlll997l : Ixitov et al.lbnod boosl ). This mapping is discon- 
tinuous at coronal current sheets (tangential discontinuities) where reconnection and en- 
ergy dissipation are particularly rapid. Regions where the mapping is extremely distorted, 
called quasi-separatrix layers, may play an equally significant role in these processes as well 



( Longcope fc Strauss! ll994J : iPriest fc Demoulinlll995l : iDemoulin et al.lll996l : iDemoulin et al. 
1997h 



While reconnection is a local process, energy storage is global, so its study requires 
a global characterization of connectivity. For example, when two active regions interact 
magnetically, there is a change in connectivity whereby new fiel d lines are forged to connect 



the p ositive polarity of one to the negative polarity of the other (ISweetlll958l : iLongcope et al. 



20051 ). One method of quantifying the global connectivity is to group photospheric footpoints 
into a number of unipolar regions. This is a tacit step in, for example, the above reference 
to the "positive polarity of the active region". Coronal field lines are then categorized by 
the regions to which their positive and negative footpoints belong. 

Global connectivity may be used to quantify the energy stored as coronal field evolves. 
Barring magnetic reconnection, emergence or submergence, the total fiux in each connection 
is preserved even as the photospheric regions move and deform. The preservation of connec- 
tion fiuxes constitutes a set o f constraints wh ich can be used to place a lower bound on the 
coronal magnetic free energy (jLongcopell200ll ). 



The process of grouping photospheri c footpoints into unipolar regions, called partition 



inq, i s natural in certain idealized models (ISweet 



2002 



1999 ) or in the fiux elements of the quiet Sun (ISchriiver et al 



1958 



Gorbachev fc Somov 19881 : Brown &: Priest 



Welsch fc Longcopd l2003l : IClose et al.l l2004l : IPeForest et al.l 120071 ) . Magnetograms of 



1997 



Hagenaai 



2001 



Parnel] 



real active regions, on the other hand, show photospheric field distributed in a complex 
pattern whose redu ction to regions is less straightforward. One algorithm de veloped by 



Barnes et al.l (120051 ) uses the vertical field in a gradient-based tessellation method (IHagenaar et al. 



19971 ). This breaks an active region into anywhere from dozens to hundreds of unipolar pho- 
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tospheric regions, depending on algorithmic parameter s. Significantly, the reg ions identified 
by the algorithm track inferred photospheric motions (ILongcope et al.l 120071 ). so their con- 
nectivity can be used to bound coronal energy. Larger numbers of regions will lead to a 
l arger number of constraints, and therefore a more restrictive lower bound on free energy 
( lLongcopell200ll : iLongcope fc Beveridgd 120071 ). 



The connectivity between photospheric regions depends entirely on the coronal magnetic 
field anchored to it. Since high spatial resolution measurements are made only at the lowest 
level of the atmosphere, such as the photosphere, it is necessary to extrapolate these data 
into the corona before connectivit y can be determiiied. T here are numerous methods for 
performing this extrapolation (see iMcClymont et al.l 119971 . for a review) and each one will 
produce a different connectivity. The most sophisticated class of methods, the non-linear 
force-free field (NLFFF) models, includes at least a half dozen variants, many of which have 



recen t ly been inter-compar ed in a series of investigations (jSchrijver et al.ll2006l : iMetcalf et al. 



20081 : ISchrijver et al. 



2008r). The conne ctivity was used as a basis of comparison in one of 



these investigation (IMetcalf et al.ll2008l ). and it was found that the different NLFF fields 
produced by these methods each induced a different connectivity. Indeed, the differences 
in connectivity tended to be large even when other metrics showed reasonable agreement 
between an extrapolation and the model field. It is therefore essential to understand how 
much the connectivity might vary under different extrapolations. This is the objective of the 
present study. 

At the opposite extreme to NLFFF models is the potential field extrapolation, which 
assumes the corona to be current-free. This assumption leads to a well-posed mathematical 
problem whose solution is relatively straightforward. Nevertheless, several versions of the 
potential field are possible depending on the treatment of the boundaries. For example, the 
magnetogram of a single active region can be extrapolated onto a finite computational grid, 
into an infinite half space, or into a spherical corona inside a source surface. Each choice has 
advantages and all are in common use. Since these fields are all different, it is to be expected 
that each will produce a different connectivity. 

In this work we will explore the difference in connectivity produced by different methods 
of field extrapolation. Since the principles of potential field extrapolation are so simple and 
well understood, we restrict our investigation to these alone. We will explore the differences 
produced by different treatments of the boundaries when making potential field extrapola- 
tions. In addition to their many other complications, sophisticated extrapolation methods, 
such as the NLFFF, must choose between these same boundary conditions. It is therefore 
worth quantifying the effect of these choices on connectivity before considering the effects of 
more complex extrapolations. 
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When a coronal field model is to be used only to compute connectivity between unipo- 
lar photospheric regions, it is possible to replace those regions by magnetic point charges. 
The result, known as a magnetic charge t opology model (MCT), is a kind of field com- 
monly used to study mag netic topology (IBaum &: Bratenahll llQSOl : [Gorbachev fc Somov 



1989 : lBrown fc Priest 
19981 



Longcope et al. 



9991) as wel l as to quantify connectivity in observed fields (ILongcope fc Silva 



20051 . 120071 ). Point magnetic charges situated on the photospheric 



level create certain unphysical artifacts, such as divergent magnetic field. These artifacts 
are absent from more traditional extrapolations which take the magnetogram itself for 
the lower boundary condition. MCT models do, however, off er the advantage that their 



topologies may be rigorously and systematically characterized ( ILongcope fc Klapperl 12002 
Beveridge &: Longcopd l2005l ). Furthermore, a potential MCT field in an unbounded half 



space takes the form of an analytic sum whose evaluation does not require a computational 
grid. 

While it is clear that point charges introduce tremendous errors in local properties of 
the magnetic field, it is not clear how much they affect its global connectivity. Away from a 
given unipolar region, the potential field will be dominated by its lowest multipole moments, 
monopole and dipole. These terms are exactly matched by a single point charge. In this work 
we quantify the difference in connectivity produced by using point charges in a potential field 
extrapolation. 

In order to make a realistic comparison of connectivities w e use a magne t ogram of an 
actual active region. This magnetogram is analyzed using the iBarnes et al.l (120051 ) parti- 
tioning algorithm with a range of different parameters. The resulting partitions consist of 
anywhere from 35 to 395 different unipolar regions. We compute connectivities between 
these regions using four different potential magnetic field extrapolations. We find that the 
connectivity of a given partition varies by no more than 15% regardless of what potential 
field extrapolations is used. Among the extrapolations, the choice of outer boundary makes 
the greatest difference. The use of point sources changes the connectivity by roughly 6% (in 
this case) regardless of how many unipolar regions are present. 

The next section reviews the process of partitioning and shows how its parameters affect 
the result. Section 3 defines the connectivity fiux and describes our method for calculating 
it. The following section describes the different extrapolation methods we explore. The 
connectivities are then compared in section 5. 
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2. Partitioning the magnetogram 



We begin with the single magnetogram from the Imaging Vector Magnetograph (IMickey et al. 



19961 . IVM) shown in Figure [H The magnetogram is of AR 8636 from 23 July 1999, and 
includes most of the flux obviously belonging to the acti ve region (AR). The i nversion of 
the spectra to produce magnetic field maps is described in lLeka &: Barnes! (120031 ). while the 
ambiguity inherent i n the observed t r ansve rse component of the field was resolved using the 
method described in ICanfield et al.l (119931 ). The three vector components of the resulting 
vector magnetic field are used to compute the vertical (i.e. radial) component in each pixel, 
Bz{x, y). These values form a 237 x 202 array of 1.1" x 1.1" pixels within the plane of the sky. 
Since the active region is relatively close to disk center, and we are using the magnetogram 
for illustration purposes, we do not project the image onto the solar surface. Instead we 
perform all analysis within the plane of the sky. 

It is evident from the cumulative histograms of positive and negative pixels, shown in 
Figure [2], that the data are dominated by positive flux. Positive pixels (B^ > 0) account 
for $+ = 5.1 X 10^ G arcsec^, while negative pixels compose less than two thirds of these 
values = 3.2 x 10^ G arcsec^. (Had the radial field been mapped to the solar surface 
the fluxes would have been 2.9 X 10^2 and i g x 10^^ Mx respectively.) It can be seen 
from the histogram that field stronger than 500 G, which accounts for ~ 2 x 10® G arcsec^, 
is much better balanced; most of the excess positive flux is weaker than this. This apparent 
imbalance probably arises from the exclusion from the IVM field of view, of an extended, 
diffuse region of negative polarity to the East. This is possibly part of an older, decaying 
AR into which 8636 emerged. There is also an excluded region of more positive polarity to 
the South of the IVM field of view. 

Some degree of flux imbalance is inevitable in any magnetogram data. Consequently, 
any method of magnetic extrapolation and any determination of magnetic connectivity, must 
somehow accommodate imbalance. Our example, with its extreme degree of imbalance, will 
bring these issues to the fore. Moreover, we show below that connections outside the AR are 
quantified more accurately and with less computation in cases of very strong imbalances. It 
is for these reasons that we select the IVM data from Figure [1] for the present study. 



The magnetogram is next subjected to a process called partitioning (IBarnes et al.ll2005 



Longcope et al.l 120071 ) whereby pixels are grouped into unipolar regions. Pixels with field 
strength below a cutoff, here set to the 3a noise level, i?th = 76.9 G, are dis carded, and the 
remaining pixels are grouped using a gradient-based tessellation scheme (IHagenaar et al. 
19971 ). This grouping uses the gradient of a field constructed by convolving with the 



-6- 



1999-07-23T1 7:38:00.0007 




-150 -100 -50 50 



Fig. 1. — The vertical field, Bz{x,y) from the IVM magnctogram of AR 8363 (grayscale). 
Axes give coordinates in arcseconds from disk center. Curves outline the regions from a 
partition with saddle level of Bgad = 100 G and smoothing oi h = 1.0 pixel. Several regions 
are labeled for future reference. 
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Fig. 2. — A cumulative histogram of the IVM magnetogram from Figure [H Sohd and dashed 
hues show the amount of magnetic flux in the plane-of-the-sky (in Gx square- arc-seconds) 
above a given field strength, in positive and negative polarity respectively. The plus and 
times mark the amount included within regions; the threshold is Sth = 76.9 G. 



kernel function 

hl2'K 

This function integrates to unity over the entire plane, and to a/1/2 within a circle of radius 
h. It therefore smoothes out fluctuations in Bz on scales smaller than ~ h. The function 
Kh was chosen because it is the Green's function for potential extrapolation upward from 
an unbounded plane to a height h. The convolution Kh * Bz therefore resembles the vertical 
field within the plane z = h. 

The gradient-based tessellation assigns a unique region label to every local maximum 
in the smoothed field \Kh * Bz\ for which \Bz\ > Bth- The label from a given maximum 
is given to all pixels which are strictly downhill with respect to \Kh * Bz\, and also have 
\Bz\ > Bth- The resulting regions are separated by areas where \Bz\ < -Bth, or by internal 
boundaries emanating from saddle points in the convolution \Kh*Bz\. The next step, called 
saddle-merging, eliminates any internal boundary at whose saddle point \Bz\ is greater than 
a value min(|i?pk|) — -Bgad, where Spk are the values at the neighboring peaks, and -Bgad is 
a parameter of the partitioning. The regions are merged by relabeling the smaller one with 
the region number of the larger. Of the remaining regions, any which have flux less than 10^ 
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G arcsec^ are discarded. 

The partition of a particular magnetogram depends critically on the parameters h and 
-Bsad; as illustrated by Figure El Increasing the smoothing- kernel width h diminishes the 
number of local maxima. The result is fewer regions which are consequently larger. Simi- 
larly, increasing -Bsad, eliminates more internal boundaries, again yielding fewer regions. The 
progression is evident in Figure [3] by scanning up the columns or rightward along the rows. 
The parameters on the upper right {h = 2.0", -Bsad = 200 G) partition the entire magne- 
togram into in = 40 regions, while those in the lower left {h = 0.1", Sgad = 10 G) partition 
it into = 395. 

Region TZa from a given partitioning is a set of pixels in which is of the same sign: 
the region is unipolar. The region may be characterized by its signed net flux 

*a = J B,{x,y)dxdy , (2) 

and its centroid location 

x'^ = J ^B,{x,y)dxdy . (3) 

Ha 

(We write integrals for mathematical clarity, but these are actually computed as sums over 
pixels in TZa multiplied by the pixel area A^^^ = 1.21 arcsec^.) 

Further characterization of a region is provided by its quadrupole moment 

.-1 



<5°j = {^i - - ) {xj - x]) B,{x, y) dx dy , (4) 



where i and j are component indices for the horizontal vectors (either 1 or 2). One measure 
of a region's horizontal extent is its radius of gyration 



9 



VQu + Q22 ■ (5) 



Its elongation may be characterized by 



where A^and are the smaller and larger eigenvalue of Q,"-. Since Qij is positive definite 
the elongation parameter will lie in the range < e < 1. An axi-symmetric flux distribution 
will have e = 0, while a very long distribution will have e ~ 1. 
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h = 0,1 



Fig. 3. — A subsection of the magnetogram (see Figure [T]) partitioned using different pa- 
rameters. The columns show the smoothing parameters h = 0.1 pixel (left), 0.5 (center) and 
2.0 (right). The rows show the saddle-merging levels i?sad = 10 G (bottom), 50 G (mid- 
dle) and 200 G (top). Note that increased smoothing removes the small positive polarity 
area surrounded by negative polarity, while increased saddle-merging simplifies it to a single 
partition without removing it. 
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For a given choice of parameters, h and -Bsad, the partitioning algorithm will break a 
magnetogram into regions with various characteristics, and r^. The partitioning may be 
characterized as a whole using a flux-weighted average of each characteristic 

|^a|/a 

a 

Figure m shows the flux- weighted averages of quantities arising from partitions with different 
parameters. Even as the parameters cover a rectangle of h and Ssad most of the averaged 
values fall on a single curve ordered by the total number of regions A^. The average elongation, 
(e) (b) decreases slightly from 0.4 to 0.3 as increases from 40 to 400. It would seem the 
flner partitionings (smaller h or smaller -Bsad) produce less elongated regions. 

The average radius of gyration (c) shows a far more pronounced decrease with increasing 
number. Clearly flner partitioning produces regions which are generally smaller. The panel 
to its right, (d), plots the distance, denoted A^, from the centroid of region a to the nearest 
neighboring centroid of either polarity. Its flux-weighted average falls along the curve ( Aq) = 
171iV~^/^. Had the centroids been scattered randomly over the ma gnetogram their nearest - 



neighbor distance would tend toward (A^) = 120A^~^/^, at large N (IKendall fc Moranlll963l ). 
The median value of A^ does, in fact, approximate this curve, but the flux- weighted mean is 
dominated by larger regions which tend to be further from their neighbors. This tendency 
leads to the larger coefficient. 



3. Connectivity 

The connectivity between regions can be found given some coronal field anchored to 
the partitioned photospheric field. There is a connection between source regions a and b if 
a coronal field line has one foot in positive region TZa and the other in negative region 
Figure [5] is a schematic depiction of all connections produced by a particular field anchored 
to the partitioning from Figure [H 

The connection between region TZa and TZb can be quantified by the connection fiux, 
ipab (the first index will always designate the positive source region). If a connection exists 
between these regions then ipab > 0; if they are unconnected then ipab = 0. All of the fiux 



^We make a distinction here between a connecti on and a domain, since it is p ossible for a pair of sources 
to be connected through more than one domain ( Beveridge &: Longcopel 20051) . Regardless of how many 
domains connect the sources, we count this as a single connection. 
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Fig. 4. — Flux-weighted averages of region characteristics for different partitions. The bot- 
tom panels plot the radius of gyration Vg (c) and the distance to nearest neighbor, A„ (d), 
both expressed in arcseconds. The top panels show the ratio rg/A„ (a) and the region elon- 
gation e (b). In each plot the value is plotted against the total number of regions in that 
partition. The different symbols denote different smoothing parameter and colors designate 
the saddle-merging level. The key for all four plots appear in the upper right of (b). An 
orange square {Bs = 100 G and h = 1.0) marks the partition from Figure [1] with its = 70 
regions. 
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Fig. 5. — A schematic diagram of the connections between source regions from Figure [H 
induced by coronal field MB. Centroids of positive and negative regions are denoted by 
pluses and times. Connections with flux in excess of 0.25% of <I>tot are represented by dark 
lines, those with less by light lines. Diamonds surround regions with connections to infinity. 
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anchored to negative region TZh must have originated in some positive region, so 



Jab , 



where S+ is the set of all positive regions. A similar expression holds for positive region a 

$a = $^^a6 , (9) 

beS- 

except that 5_ includes an extra source to account for field lines extending to "infinity" 
due to the net positive imbalance. This new source has flux |$oo| = '^'+ — I'^'-l, in order to 
account for all flux which cannot close in a photospheric negative source. With the inclusion 
of infinity as a negative source, sources of each sign will have the same net flux 

$tot = Y.^'^ = -^2^"= E E ^-^^ • (10) 

aes+ beS- ae<s+ beS- 

This is the total amount of flux in the field. 



The connection flux ipab can be estimated using a Monte Carlo method (IBarnes et al. 



20051 ). A number of footpoints are selected randomly within positive source region TZa- 
Each field line is then followed to its other end (or to a distance from which which it will 
not return). The number which end at region TZb is designated rriab- Field lines are then 
randomly initialized from negative regions (except infinity) and traced backward to their 
positive source region. The number initiated in TZfy whic h "terminate" in T Zn are denoted 
rriba- The Bayesian estimate of the connection flux ipab (IBarnes et al.l l2005l ) combines the 
information from tracing in both directions as 

, _ mab + rriba , 
- nj<l>^ + nbm ■ ^ ' 

Since no points were initiated at infinity, Ub = and rriba = for that case {b = oo). The 
estimate then reduces to V'afe = {^ab/na)^a- 

Expression (fTTl) is a Monte Carlo estimate of the actual connection flux so it will in- 
clude some statistical error. This error can be estimated using the (approximately) Poisson 
statistics of the counts ruab and rriba- If an actual connection includes very little flux it is 
possible that none of the randomly generated field lines will sample it and rriab = rriba = 0. 
The probability that an actual connection will be erroneously missed for this reason is 

Pmiss(V^a6) = e-'^-/^" , (12) 

where (po = l/iria/^a + nb/\^b\)- Naturally the use of more points (i.e. larger values of ria 
and rib) will make this increasingly unlikely. Nevertheless, there will always be a possibility 
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that some number of very small connections have remained undetected. For one particular 
case (later defined as field PH) a topological analysis of the connections revealed that 10 of 
109 connections were erroneously missed by the Monte Carlo method. 

4. Coronal fields 

4.1. Different potential fields 

In order to define connectivity it is necessary to compute the entire coronal field an- 
chored to the photospheric field Bz. Although we have chosen to restrict consideration to 
potential fields, V x B = 0, there are several different ways to compute a potential field 
from magnetogram data. We consider a variety of these fields and study the effect on the 
connectivity of the differing fields. 

One option is to compute the potential field within a rectangular x Ly x Lz, box, 
B, with a lower boundary at the magnetogram, z = 0. The four lateral boundaries are 
perfect conductors (n • B = 0) positioned along the edges of the magnetogram. The vertical 
field at the lower boundary, Bz{x, y, 0), is taken from the magnetogram and is therefore not 
balanced. An equal net flux must cross the upper boundary or no solution would be possible 
for which V ■ B = 0. We achieve this with a uniform vertical fleld along the upper boundary: 
Bz{x, y, Lz) = \^oc\/ {L^Ly). A field line crossing this upper boundary is designated as a 
connection to infinity. We choose to place this upper boundary at — 220", approximately 
equal to Ly, and slightly less than L^.. 

The alternative to the computational box B is to use a coronal field extending throughout 
the entire half-space z > 0. Such an unbounded field is computed, in principle, by convolving 
the field at z = with a Green's function for a point magnetic charge aX z = 0. In practice 
we compute either a portion of the field on a grid or compute it along a field line as we trace 
it. We distinguish between the box boundary and the half-space using superscripts B and 
H respectively. 

For the purposes of computing connections between unipolar regions it is possible to 
replace each region with a point source. Region TZa is replaced with a magnetic charge of 
strength = $a/27r at position x" on the photospheric plane, z = 0. This matches a 
multipolc approximation of the potential field from TZa up to the dipole term, thus it is 
expected to be accurate at distances |x — x"| » r^. Fields computed using point sources 
will be assigned a superscript P. With this simplification the convolution required for the 
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half-space computation becomes 

a 

where the sum is over all N sources, not including infinity. 

An alternative to photospheric point sources is to compute a potential field matching 
the magnetogram pixel-for-pixel. Field vectors are computed on a three-dimensional uniform 
cartesian grid, for example within B. In order that each field line have a defined connectivity 
magnetogram pixels which do not belong to a source region are set to zero. The resulting 
magnetograms differ slightly for different partitionings, but for most cases = 4.7 x 10^ 
G arcsec^ and $_ = 2.7 x 10^ G arcsec^. Fields anchored in this way are designated by a 
superscript M. 

The field B^^^^^^ bounded by conducting boundaries, is readily computed on the Carte- 
sian grid using Fourier methods. The extrapolation into the half-space, B^^-^^ is done onto 
a larger Cartesian grid which we call Qmh- The field at any grid point can be found by 
convolving the entire magnetogram with the half-space Green's function. Performing the 
convolution for every grid-point is very time-consuming so we use it only for points along 
the boundary OQmh- We then use the efficient Fourier techniques to compute the inte- 
rior potential field matching these boundary values. It is noteworthy that fiux crosses the 
boundaries dQMH, so these are not genuine boundaries of the field. 



4.2. Calculating the connection fluxes 

For each type of coronal field, connection fiuxes ipab are computed using the Monte Carlo 
estimate in eq. (ITT]) . Approximately = |$a|/(20o) field lines are initiated on source a, 
where we have chosen (po = 10 G arcsec^. The smallest fiux reported will be ipab — (po, for 
example when mab = 1 and rriba = 0. With our choice of (po, every source will have at least 
50 lines, and the largest will have npoi — 4 x 10^. One complete estimate, such as the one 
shown in Figure O requires ($+ + |$_|)/20 ~ 3.7 x 10^ field lines be traced. 

Field line initialization is different for the cases anchored to point sources than for those 
anchored to magnetograms. In the point-source cases, the Ua points for a given source are 
randomly generated with a uniform distribution over a very small hemisphere centered at 
the point source. The radius of the hemisphere is set to be small enough that the magnetic 
field is directed roughly radially outward (inward) from the positive (negative) source. 

For magnetogram cases the extended region TZa is a set of Pa pixels. First, a list 
of Ha pixels are randomly generated so as to sample each pixel in TZa with probability 
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proportional to its field strength (p = ylpix-Bz/$a)- Since Ha ^ -Pa, typically, the list will 
include many duplicate pixels. For each random pixel (including all duplications) an initial 
point is randomly generated with a uniform distribution over the pixel. Thus no two field 
lines from the same pixel will begin at the same point. The subsequent field line integration 
uses tri- linear interpolation to calculate B(x), so different points within the same pixel will 
belong to different field lines. 

Using these methods we estimate the connectivity of a given field according to Equation 
f llip . As an illustration, consider the field B*^^^^) from the partition shown in Figure [H Our 
estimate, shown in Figure O includes 198 different connections between its = 70 sources 
(71 including oo). These connections are quantified by their connection fiuxes, ipab, plotted 
in Figure O The connections to a given source fall along a vertical line below the diagonal, 
V'afe = (dotted) in the lower panel. The number of connections to that source, called its 
degree da, is plotted above its fiux in the upper panel. 

Some sources, such as POl or A^19, have only a single connection {da = 1) and are called 
leaves. In Figure leaves appear at the end of a single line (A^19) or as an isolated diamond 
(POl). In the lower panel of Figure El the fiuxes of a leaf connection naturally fall on the 
dotted diagonal since all the fiux from that source belongs to that single connection. 

In contrast to the leaf connections there are several sources, such as A^Ol with many 
connections (rfjvoi = 27). Given C = 198 connections to + 1 = 71 sources, the average 
source must connect to da = 2C/{N + 1) = 5.65 sources: the value marked by the dotted line 
in the upper panel. (The factor of two arises from the fact that each connection is incident on 
two different sources: one positive and one negative). There is a notable tendency for larger 
sources, especially larger negative sources, to have more connections. As a consequence of 
this tendency, the fiux-weighted average. Equation ([7]), of the degree is (da) = 9.35 in this 
case. 



4.3. Connections to infinity 

In each different magnetic field there is open fiux, represented by field lines connected 
to infinity (formally a negative source). In the MB and PB fields, an open field line is one 
that terminates at the upper surface of the box, z = L^. The 34 positive sources enclosed 
by diamonds in Figure O are connected to infinity in the MB field for that partition. These 
form the connections whose fiuxes fall in the right-most vertical row of x s in Figure [61 

The PH field occupies the entire half space and open fiux truly extends outward indef- 
initely. Far from the AR the field resembles that from a single point charge I^ooI/Stt; field 
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Fig. 6. — The connection fluxes ipab from the field depicted in Figure O The bottom panel 
plots the fiux ipab against the fiux of its sources, (as a +) and \^b\ (as a x). The sloped 
and horizontal dotted lines show the maximum and minimum possible values: ipab = and 
00 respectively. Colors are used to differentiate the different sources, whose connections fall 
along a vertical line. On each axis the fluxes are normalized to the total, $tot- The top panel 
plots the number of connections to each source, da- Several of the sources are labeled for 
reference. The horizontal dotted line shows the mean degree, da = 5.65. 
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lines go outward approximately radially. There is a single separatrix surface dividing closed 
from open flux, and once a field line has been integrated far enough to establish that it lies 
outside this surface, it may be designated as a connection to infinity. 

The separatrix between open and closed field is a dome anchored to a null point located 
~ 327 arcsec from the center of the AR (the triangle in Figure [7]). The footprint of the dome 
passes along a series of spines (solid curves) linking positive sources. These sources connect 
to infinity as well as at least one other negative source, A^02, inside the dome. Positive 
sources outside the dome link only to infinity. 




Fig. 7. — The footprint of the separatrix dividing closed from open flux (dashed and solid 
curves) in the field B'^'^^^ from the partition of Figure [H Positive and negative sources are 
denoted by +s and xs respectively. Triangles are magnetic null points. The triangle at the 
far left, lying along the dashed curve, is the negative null whose fan surface is the separatrix. 
The center of unsigned flux is denoted by a diamond, and the dipole moment about this 
point is indicated by an arrow. The broken rectangle shows the extent of the magnetogram, 
and the larger dotted one is the bottom surface of Qmh- 

The dipole moment, computed about the center of unsigned flux (diamond in Figure [7]), 
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is fi = 5.1 X 10'^ G arcsec^, directed 7° below x as shown in the figure. A far field with this 
dipole and the net charge, goo = +|'^'oo|/27r will vanish at one point located a distance 

ro = 2— ~ 318arcsec , (14) 
goo 

from the center of unsigned charge. For a region with positive net fiux, as we have, the null 
is situated in the direction opposite to the dipole moment. Clearly this is the approximate 
location of the null point whose separatrix divides open from closed fiux in the PH field. 
Had the region been more balanced, goo would have been smaller and the separation of open 
from closed fiux would have occurred much farther out. 

There are 22 positive sources linked to infinity by dint of lying on or outside the sepa- 
ratrix in field PH. In contrast, the MB, field has 34 sources connected to oo, even though 
there is the same amount of open fiux, |$oo|; in each case. Comparing the labeled sources 
in Figures [5] and U\ shows that different sources are so connected in the MB and the PH 
fields. These connections are just some of the differences between the two fields, explored 
further in the next section. 

Establishing and quantifying connections to infinity is particularly challenging for the 
MH field. Although this field formally extends throughout the half-space, it is known only 
on a Cartesian grid covering Qmh- A field line can therefore be followed to the boundary of 
Qmh but no farther. It is, in principle, possible for field at dQuH to be directed both inward 
and outward. When this is the case there must be some field lines which leave the volume 
where i?„ > and return where Bn < 0. Naturally these field lines cannot be followed, so 
their fiux cannot be correctly assigned to a connection. Due to the large fiux imbalance in 
our magnetogram we were able to select a volume Qmh for which Bn > on all outer surfaces 
(see Figure [7]). Thus any field line encountering the boundary is necessarily connected to 
infinity. 

In order to assure Bn > on the outer boundary it is necessary (but not sufficient) 
that Qmh enclose the separatrix dome in the field 'B^'^^K This surface will closely resemble 
(but not exactly match) the separatrix dome from B*^'^^-' since both fields approach the 
same far-field form. It is evident from Figure [7] that the base of Qmh does enclose the 
latter separatrix. This requirement means Qmh must be considerably larger than B, so 
computations in B*^*^'^^ are much more expensive than for B*^*^^) or B^^^\ If the active 
region had had better fiux balance then Qmh would need to be still larger and computation 
would have become prohibitively expensive. 
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5. Comparisons 

5.1. Different fields from a single partition 

Four different coronal fields can be generated from a single partition in the fashion 
described above. Each field will contain the same total flux, $tot; interconnecting the same 
+ 1 sources. The connections will not, however, be the same for the different fields. Figure 
[H] shows the connection fluxes ipab^^ (the ones from Figure E]) plotted against those induced 
by B^^^\ There are C^^'^^^ = 198 connections in the former and only C^^^^ = 184 in the 
latter. Moreover, there are 48 connections in MB which do not appear in PB; these appear 
as diamonds along a vertical line in the central panel. Similarly 34 connections in PB not 
present in MB form the horizontal row of diamonds. We refer to either of these as singlet 
connections. The remaining 150 connections, common to both fields, are plotted as +s in 
the central panel. 

The tendency for common connections (+s) to cluster about the diagonal, especially at 
the upper right, shows that connections have similar fluxes in both fields. The connections 
falling outside the dotted diagonals (i.e. disagreeing by more than a factor of two) are over- 
whelmingly dominated by smaller connections: ipab < 10^^$tot- These small connections 
also compose almost all of the singlet connections in either field. Indeed, a great many of the 
singlets are so small {ipat — 0o) that they had a significant probability of going undetected 
even in the field where they were found. These tiny connections account for most of the 
spikes at the small-flux end of each histogram. 

The statistical errors from the Monte Carlo calculations are relatively large for small 
fluxes (found by only a few field lines). On the logarithmic plot, like Figured the error bars 
are largest at the bottom or left. Above a value of V'afe — 10~^$tot statistical errors are less 
that 0.1% and errors bars are smaller than the symbols. 

The impression given by a comparison such as Figure [8] is that in spite of their differences 
the two coronal fields induce connections which are largely in agreement. The differences 
appear mostly in the very small connections. While these are small, most of them lie well 
above the detection limit, 0o; and thus represent genuine differences. The fact that the 
differences are in small connections suggests that they will not be of great importance to a 
model of the field. 

In order to weight the most significant flux differences we focus on the difference 

A'^u=i.r-i'T'' ■ (15) 

Connections with positive difference {Aipab > 0) are those for which MB has excess flux 
relative to PB; these appear above the diagonal in Figure [H Using Equation ffTOl) we can 
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Fig. 8. — The fluxes of identical connections from two different flelds, MB (vertical axis) 
and PB (horizontal axis) plotted against one another. Fluxes are normalized to the total 
^tot- Vertical and horizontal dotted lines show the minimum flux, (pQ. Singlet connections 
(absent from one of the fields) appear as diamonds below or to the left of these limits. 
For example, the vertical line of diamonds are absent from PB but present in MB with 
fluxes indicated by their position. The diagonal dotted lines mark i^ab^^ ~ ^'^if^^ (upper) 
and ip'^J^^'' = /2 (lower). Statistical errors are indicated on a few representative points. 
Plotted along the right and top are histograms of In ipah for that fleld. The number appearing 
in the histogram gives the total number of connections. 
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show that 

so there will be as much flux in connections with excess {^il^ab > 0) as in connections with 
deficit {Aipab < 0). Figure shows cumulative histograms of the flux differences of each sign. 
In each case the total discrepancy is 6.80% of $tot- The connections are sorted by decreasing 
magnitude, so the histograms rise sharply at flrst. Ten to twelve connections account for 
half the total discrepancy of each sign (shown by diamonds). The rest of the discrepancy 
occurs in the hundred or so other connections. 
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Fig. 9. — Cumulative histograms of the flux differences Atpab- Those with Aipab > are 
graphed on the right and the ones with Atpab < are on the left. The differences are sorted 
by decreasing magnitude, so the upper curves plot the total difference for the n largest 
differences of that sign. The right axis gives this flux in units of G arcsec^ while the left 
axis gives it in units of $tot- The lower curves are for only the singlet connections. A + 
and X appears at the extreme values of the total and singlet curves respectively. These 
extreme values are also projected to the ordinate and abscissa with the same symbol. The 
sole exception is that the x on the left gives the sum of the total singlet fluxes of both signs, 
while the curves show the sums for each sign separately. A o shows half the total discrepancy 
for each sign of Aipab- 
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For a singlet connection, one of the terms on the right of Equation f|T5|) will vanish. The 
magnitude of the difference will therefore equal the other term. The lower curves in Figure [9] 
show histograms formed from the singlet connections alone. The right curve accumulates the 
differences in the 48 singlet connections in MB, while the left accumulates the 34 singlets 
in PB (for which Aipab < 0). Combining their totals yields 1.05%$tot, so singlets contribute 
only a small fraction to the overall flux discrepancy (6.80%). 

The inclusion of two other kinds of coronal field leads to six different pairwise compar- 
isons of the kind just used. Figure [TD] shows cumulative histograms like those in Figure [HI for 
all six possible flux differences. All histograms have the same T shape as the ones in Figure 
m a few connections account for the majority of the discrepancy. Variation among the three 
fields seems to be mostly a matter of degree. 




Fig. 10. — A superposition of cumulative histograms comparing all four different kinds of 
fields. Each comparison is the same as Figure [H but shown in a different color. The color 
code appears along the top. 

The largest discrepancies are between fields differing both in the their anchoring (mag- 
netograms versus points) and their boundaries (box versus half-space): for example between 
PB and MH. The fields differing in only one respect, either anchoring or boundary, have his- 
tograms in a cluster below the other two. This tendency is true for both the total histograms 



-24- 



(upper curves) and singlet histograms (lower curves). 

Making the same plots for many other partitions we find no exceptions to this T shape. 
The histograms therefore differ mostly in their magnitude, which may be summarized by the 
maximum (marked by a + along the left axes in Figures [9] and [TOj) . We denote this value 

aG5+ 6g5_ aG5+ 6g5_ 

for a comparison between B*^"^) and B^^^ for the same partition. (Here X and Y stand for any 
of MB, PB, MH or PH.) Thus Figure[Tlis summarized by the values /\^pb-mh = 10.8%, 
A-^MB-PB = 8.6%, A-^MB^PB = 6.8%, and so forth (in units of $tot)- 

(X) 

All the connection fiuxes are calculated by Monte Carlo methods, so the values of t/'^^ 
and 'i/'^P include statistical errors. The sum in ( |T71) will therefore be biased upward, and 
even an estimate of A'^x-x will be positive provided it uses two different estimates of the 
connectivities from B*^"^'. (To see this note that the sum in Equation [l7j is over numbers 
which are never negative and are usually positive due to errors.) We can subtract the 
expected bias, assuming errors in Aipab to have Gaussian distributions, following a procedure 
described in an appendix. Doing so for A'^mb-mb, for example, yields a number consistent 
with zero. Doing so for the red curve in Figure [TO] yields A'^mb-pb = 6.6% ± 0.25% (the 
value on the curve, A^^^mb-pb = 6.8%, is therefore biased upward by ~ 0.2%). Thus there 
is a true difference between field MB and PB, beyond that caused by statistical errors. 

The quantity A\E'x-y can be considered the connectivity distance between the coronal 
field models B*^"^) and B*-^). When A'^x-y = the fields are identical, at least with respect 
to their connectivities. For a given partition there are four different fields separated by six 
distances. The fields can be represented as vertices of a tetrahedron in a three-dimensional 
space. The left column of Figure [TT] shows two views of the tetrahedron formed from the 
histograms in Figure [TOl The distances used in these plots have biases removed. 

If all six distances were exactly the same, they would form a regular tetrahedron. In fact 
two of the distances, A'^pb-mh and A"^ mb-ph, are largest, leading to the extremely fiat 
tetrahedron shown in Figure [Til The fiattened shape is approximately a quadrilateral lying 
in a plane, with the two large distances forming its diagonals. The sides of the quadrilateral 
separate vertices (fields) differing in only one respect. The face-on views of the quadrilateral 
(bottom) are oriented so that horizontal edges separate fields of different anchoring (M vs. 
P) while vertical edges separate fields of different boundaries (B vs. H). 

Performing the same analysis for different partitions gives distances with similar prop- 
erties, as exemplified by the center and right tetrahedra in Figure [TTJ Fields differing in 
both anchoring and boundaries are furthest apart, so the tetrahedron is flattened into a 
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Fig. 11. — Tetrahedra formed from the connectivity distances between four fields anchored to 
different partitions. Each tetrahedron is oriented so that MB, PB and MH he in a plane. 
The view perpendicular to that plane is shown below and tangent to the plane is shown 
above, with the same scale. The tetrahedron on the left is from the partition of Figure [H 
whose distances appear in Figure [TUl Distances are in percentages of <l>tot and horizontal and 
vertical reference bars are shown (5% and 10% respectively). The center and right columns 
are tetrahedra for partitions with N = 108 and N = 395 regions. The vertices, labeled on 
the left, are in the same orientation for all three cases. 
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quadrilateral. Furthermore both instances of a particular difference, such as the anchoring 
(the horizontal sides in the lower view), result in a similar distance. This makes the quadri- 
lateral into an approximate parallelepiped. For the cases with more sources (center and 
right) changing boundaries makes the largest differences, so the parallelepipeds tend to be 
taller than they are wide. Finally, the similarity of the diagonals pushes the shape toward a 
rectangle. 

5.2. Variation in partitions 

To see how the above comparisons are affected by different partitioning parameters 
we perform Monte Carlo calculations for fields from different partitions. The smoothing 
parameter h is varied from 0.1 to 2.0 arcsecs, and i?sad from 10 G to 500 G. 

Once again we find that the different partitioning can be approximately ordered by 
the number of regions, N . Figure [T^] shows the average degree of a source region in each 
of the fields over a range of partitioning. The lower points show the un-weighted average, 
da = 2C/{N+1), as illustrated in FigureO This quantity is very similar, da — 6, for all fields 
and all levels of partitioning. It seems that the number of connections scales as C ~ 3N 
even as the number of possible connections goes as ~ A^^. In contrast, the flux weighted 
average, (da), does appear to increase with the number of regions, although perhaps at a 
power less than ~ A^. This shows that the largest sources connect to more sources as they 
become available. 

Comparisons from the previous section, between all four kinds of field, suggested that 
it is sufficient to compare only three. The six possible comparisons between all four fields 
were visualized as a tetrahedron of distances. It was found, however, that these tended to 
form a flat rectangle, well characterized by two of its sides. Taking advantage of this we 
consider the fields MB, PB and PH for a large number of different partitions. Among the 
three comparisons there is one differing only by anchoring [MB versus PB), one differing 
only by boundaries {PB versus PH) and one differing in both respects [MB versus PH). 

Finer partitioning (i.e. smaller values of h or -Bsad) result in more source regions, A^. 
According to the lower curve in Figure [12] these sources are interconnected in a propor- 
tionately large number of ways, C, regardless of which coronal field is used. One expects 
that subdividing the same total flux, <l>tot, into a larger number of pieces would yield bigger 
discrepancies, A^'x-y- Figure [13] shows that this expectation is borne out when comparing 
fields with different boundaries, [PB to PH, +s or MB to PH, *s). These curves trend 
upward as A^ increases. 
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Fig. 12. — The average degree (number of connections) per source for various fields. Different 
partitions are plotted versus the number of sources, as symbols. The average of all fields 
is plotted as a solid curve. The bottom set are the average, da = 2C/{N + 1), while the 
upper are the flux weighted average, (da) from Equation ([7]). 
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Fig. 13. — The total connectivity difference, A'^x-y, defined in Equation (fTTI) for dif- 
ferent fields and different partitions. Comparisons are between fields with different an- 
choring {MB — PB) X, different boundaries {PH — PB) +, and differing in both 
[MB — PH) *. All differences are plotted as a percentage of $tot- The dashed line shows 
(A\l/|^5_p^ + A\l'pj:^_p^)^/^ for comparison. 
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In keeping with the results of the previous sub-section, the cases which differ in both re- 
spects (*) are separated by the greatest distance. The quantity (A^|^^_p^ + A^p^_p^)^/^, 
plotted as a dashed line, appears to match the the asterisks well. This fit corroborates the 
observation from the previous section that the connectivity distances formed a flat rectangle; 
the dashed curve is the hypotenuse of a right triangle formed from the other two distances. 

A truly remarkable feature of Figure [13] is that the two bounded fields, differing only 
in their photospheric anchoring, MB and PB, differ by a relatively small amount (~ 6%) 
which does not change with finer partitioning. This corroborates the tendency observed in 
Figure [11] for the rectangles to grow taller with increasing N, without growing wider. So 
while there are ever more connections being compared, and more positive terms in Equation 
f [T7|) . the total difference does not seem to change. It seems that the difference between using 
point sources or using the actual magnetogram, is about 6% of the connectivity. 

5.3. Variations in box size 

We can further explore the effect of outer boundaries by increasing the size of the 
conducting box. To do this the box B is augmented by layers of equal width, w, along 
all boundaries except the bottom {z = 0). This new domain, called B~^, has conducting 
boundaries on the four lateral walls and a uniform field at the top boundary z = + w. 
For the field, PB^, the lower boundary has the same point sources located at the same 
positions within the central x Ly square. The field is computed on a grid with cubic 
pixels, 1.1" on side, just as in the PB field. 

Connectivities are computed in fields PB^ anchored to point sources from the partition 
with = 174 sources (-Bsat = 100 and h = 0.1"). This is computed for different boundary 
layer widths, w, and the results are compared to the PB and PH fields. For vanishing layer 
width {w = 0) the "augmented" volume corresponds to B so A'ifpB-pB+ = 0. In the other 
limit, w — > oo, Pi?+ — > PH and A'^ph-pb+ 0. Figure [TJ] shows the continuous transition 
between these limits. For the case with a border width w ~ 70" the field PB~^ has become 
equally dissimilar to both other fields, PB and PH. Borders wider than this yield a field 
still closer to that of the half-space (PH) . 

6. Discussion 

Connectivity characterizes coronal magnetic field in a manner useful for understanding 
energy release and reconnection. It is possible to quantify the connectivity of an active region 
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Fig. 14. — Comparison of fields PB^ to PB (dashed) and PH (solid) for the partition with 
= 174 regions. The symbols above the line are the raw differences, !\'^ pb-pb+ (+) and 
A^PH-PB+ (><)• The corresponding hnes are corrected for statistical bias. The dotted line 
shows the value I^^pb-mb for comparison. 

based on a single photospheric magnetogram. It is necessary to first construct the coronal 
magnetic field using some kind of extrapolation and then to partition the magnetogram into 
unipolar regions. Techniques for accomplishing each of these steps have been developed 
and are in common use. The foregoing work has presented new techniques for quantifying 
the differences in connectivity for different fields anchored to the same set of sources. This 
comparison was used to assess which steps the computed connectivity is most sensitive to. 

Our comparisons show that the connectivity is relatively insensitive to variations in the 
methods of extrapolation or photospheric anchoring. Among the cases we considered, the 

greatest discrepancy between any two fields was 15% of the total fiux. That is to say that one 
field may be converted into the other, at least in terms of connectivity, by reconnecting 15% 
of its field lines. The majority of the difference occurred in a small number of connections 
present in both fields but with different fluxes. The vast majority of connections were 
common to both flelds, however, there were instances (singletons) of connections present 
in one field and not in the other. These topological differences were found to occur most 
frequently in small connections which taken together accounted for a small part of the overall 
flux difference. 
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It is not immediately clear how large a difference one could expect between any two ar- 
bitrary fields anchored to the same set of photospheric regions. Equations ([8]) and Qj place 
numerous constraints on the possible connectivities, which could render 100% difference im- 
possible. It is worth considering a few artificial connectivities for the purpose of comparison. 
One class of connectivities are those minimizing or maximizing the informational entropy 
function, 



subject to the constraints from Equations ([8]) and iQ. The partition from Figure [T] has 
46 positive regions and 25 negative regions, including oo. The informational entropy is 
maximized [H = 4.97) by connecting the sources in all 46 x 25 = 1150 possible ways 
(V'afe = ^a|$fe|/^tot)- It is minimized {H = 3.07) by a set of 66 connections. These extreme 
cases differ from one another by A\E' = 76%. All versions of the potential field extrapolations 
have very similar entropies, H ~ 3.82 and differ from the minimum and maximum entropy 
connectivities by similar amounts: A\E' ^ 62% and A\E' ~ 58% respectively. It seems that 
potential field extrapolations are far more similar to one another then to these particular 
fields. 

It would be better to compare to connectivities generated in a more realistic fashion 
than to those extremizing an ad hoc function. We could ask, for example, how different 
is the connectivity of a potential field from that of a NLFFF extrapolated from the vector 
magnetogram; or we could compare the potential field from one time to that at another time 
(provided the magnetograms are partitioned into equivalent regions). Comparisons of this 
kind promise insight into energetics and reconnection in real coronal fields, and will be the 
topic of future investigation. In order to gain this physical insight, however, it is essential to 
know the level of difference that arises from non-physical variations such as in anchoring or 
boundary conditions alone. The present study provides that important point of reference, 
and will therefore serve as a baseline in the future studies. 

Connectivity difference has proven itself a useful metric for quantifying discrepancy 
between different coronal extrapolations from the same data. The presence or absence of 
conducting boundaries are found to have the greatest effect on the connectivities of a poten- 
tial field. Figure corroborates the expectation that more distant boundaries give a better 
approximation of no boundaries at all. The 260" x 221" magnetogram considered here can 
be expanded to four times the area, by padding with w = 120" on all sides, to produce a field 
only 4% different from that in an infinite half-space. It is possible that proportionately more 
padding would be required for magnetograms with better flux balance, since these would 
have longer closed loops. Expression f|T^ . giving the extent of the closed field, is inversely 
proportional to the degree of balance. A future study will seek a general expression for the 
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required padding to a given magnetogram. 

Alternatively, the field from an infinite half-space field, Mif, can be computed on a 
grid after using a Green's function to compute it on all lateral boundaries. Unless this grid 
is large enough, ~ tq from Equation (IT^ . there will be connections extending outside the 
grid which cannot be followed precisely. It might be possible to follow them approximately 
with a gridless field, such as PH, in a kind of hybrid method. Alternatively, it appears that 
the point source approximation alone, PH, is a relatively accurate approximation to MH 
(differing by roughly 5.5% in our case) for which a grid is not necessary. 

We find that only a small error (5%-6%) is incurred when connectivity is computed 
using a simplified, gridless extrapolation from magnetic point charges [PH or PB) in place 
of more traditional extrapolation from a full magnetogram [MB or MH). These point- 
charge models differ significantly from the actual field: for example they are singular at the 
charges. The connectivity, however, seems only mildly sensitive to these local differences. 
Moreover the connectivity difference does not increase even as the number of source regions, 
and therefore the number of connections, increases. This seems explicable by the fact that 
the point source anchoring differs from the magnetogram only within a small neighborhood 
of the charge. The differences may therefore be confined to a layer z < (r^) which shrinks 
with finer partitioning. 

Large scale connectivity is defined in terms of unipolar source regions into which the 
photospheric field (magnetogram) is partitioned. Variation of parameters controlling this 
partitioning leads to significant changes in the source regions and therefore the connectivity. 
At least for the two parameters whose variation we explored, h and -Bsad, most differences 
could be ordered just by the number of regions A^. The sizes, shapes and interrelation of 
regions appears to scale with A^, as did most properties of the potential field connectivity. 

Smaller values of partitioning parameters h or Bgad led to finer partitioning, with more 
sources and therefore more connections. Remarkably we found that the total number of 
connections increased only as A^ rather than as A^^ like the number of possible connections. 
Indeed, we found that this particular active region had approximately 6 connections to every 
source independent of partitioning parameters. Further study will reveal whether this trend 
persists in other active regions. 

The work was supported by a grant from NASA's Living with a Star TR&T program. 
GB was partially supported for this work by the Air Force Office of Scientific Research under 
contract FA9550-06-C-0019. 
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A. Estimating and correcting bias in the absolute value 

Consider an unknown quantity x whose measurement, x, includes an additive Gaussian 
error of known variance o"^. The absolute value of the measurement, \x\, is an estimate of 
whose expectation is 

= J \x + e\p{e)de = \x\ + \J^(Te-''' - \x\eiic{\x\/ V2a) , (Al) 

where erfc is the complementary error function. The second and third terms on the right 
represent a bias in the estimate of 

Eb{x) = y|ae-"'/'"' - \x\eTfc{\x\/V2a) , (A2) 

since its expectation does not vanish. For values |x| ^ a the bias error is extremely small 
(~ cre"^'^/^""^), and {\x\) ~ For small magnitudes, on the other hand <^ a), the 
estimate will be dominated by the magnitude of the Gaussian noise so ~ 0.8 a. 

The actual bias, Eb, depends on the quantity |x| whose value we are trying to learn from 
the measurement x. We cannot, therefore, subtract the exact bias from the measurement. We 
must instead construct a function of the measured value, x, whose expectation approximates 
Eh{x). This function has a discontinuous derivative at x = 0, due to its second term, 
and is therefore very difficult to reproduce in the expectation of a function of x. The 
expectation of a given function f{x) can be expressed as the convolution of / with the 
Gaussian distribution of e. This convolution effectively blurs f{x) over a scale a, thereby 
smoothing out discontinuities. 

Because of the blurring property described above subtracting Eb{x) would remove a 
broader function from the expectation of the estimate. We seek instead a function more 
sharply peaked, whose convolution will be limited to |x| < a. Following this logic we propose 
the function 

B{x,a) ^ ^e-("^MV2 ^ (A3) 
V 27r 

where a is an adjustable parameter defining the width. The expectation of the function 
resembles the first term in (1A2I) and has the same integral as the actual error 



{B)dx = I Eb{x)dx = , (A5) 

oo 
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independent of a. 

The bias estimator in Equation flA3p is limited to \x\ < a/a, and in the hmit a — > cxd 
it becomes a Dirac 5-function: B{x) —>■ a'^6{x). It is natural that in the (5-function limit the 
expectation, (lA4p . is simply the distribution of noise. For large a the bias correction will 
only rarely be non-negligible; then it will be large to compensate for the numerous times it 
was negligible. This will introduce additional variance to inferred values. Adopting instead 
a small value of a will subtract a small amount from more measurements, owing to the much 
broader scope of (B). We have found a = 3 to be a reasonable all-around compromise since 
its scope is very narrow while introducing little additional variance. 



Fig. 15. — Effects of the bias estimator B{x) on estimates of |a;|. (a) The bias error as 
a function of |a;| (in units of a). Diamonds show Eii{x), the bias in Dashed lines 

show Ef,{B) for choices a = 1/3 and a = 1. The pluses are for the value used in practice: 
a = 3. Symbols are the results of 10^ trials with Gaussian noise, while the lines are the 
analytical functions from Equations ( lA2p and ( ]A3I) . (b) The computed values with and 
without compensation. Diamonds and plusses are, as in (a), the results of averaging |x| and 
|x| — B{x) respectively. Error bars are the standard deviation in the 10^ trials. The solid 
line shows the actual value |x| — B{x) versus \x\. 

The problem we face is to compute a sum of magnitudes of measurements, Xi, of different 
underlying values, Xi. We estimate this by the sum 

(A6) 

i i 

where B{x, a) is defined by Equation (1A3I) with a = 3. Terms of the sum where \xi\ < 0.46crj 
are negative in order to remove the bias. These negative terms, as well as positive terms 
where |5;| < 0.7a, actually underestimate the bias on average (see Figure [T5l) . In order 
to compensate, those values in the range 0.7cr < < 2.5a over-estimate it on average. 
Provided the underlying values, distributed relatively uniformly within the range 

|a;| < 2.5(7, the underestimates and overestimates will balance one another, due to Equation 



~ '^\\xi\-B{xi,ai) 
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f lA5l) . thereby removing the bias precisely. Even when this is not the case, the bias error at 
an individual value of x is reduced by at least half. 
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